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We present a theoretical analysis of the evaporative cooling of an atomic beam prop- 
agating in a magnetic guide. Cooling is provided by transverse evaporation. The atomic 
dynamics inside the guide is analyzed by solving the Boltzmann equation with two different 
approaches: an approximate analytical ansatz and a Monte-Carlo simulation. Within their 
domain of validity, these two methods are found to be in very good agreement with each 
other. They allow us to determine how the phase-space density and the flux of the beam 
vary along its direction of propagation. We find a significant increase for the phase-space 
density along the guide for realistic experimental parameters. By extrapolation, we estimate 
the length of the beam needed to reach quantum degeneracy. 

I. INTRODUCTION 

Forced evaporative cooling of trapped gases is a very powerful technique to increase the phase-space density 
of an ensemble of atoms up to quantum degeneracy |0] . Particles with a sufficiently large energy (typically 
5 times the thermal energy fcsT) are eliminated. Elastic collisions between the remaining particles restore 
thermal equilibrium with a lower temperature and an increased phase-space density. The most prominent 
success of evaporative cooling is the achievement of Bose-Einstein condensation of alkali and hydrogen atomic 
gases 

In this paper we present another possible application of evaporative cooling. We consider atoms moving 
freely along the 2;-axis and transversely confined by a magnetic field gradient, which provides a harmonic 
potential in the x—y plane (see fig. [^). The atoms are injected in the plane z = with given flux, longitudinal 
and transverse velocity and spatial distributions. We suppose that a radio-frequency fleld is applied so that 
only atoms inside a tube-shaped region remain trapped, those outside are evaporated. The bias magnetic 
field along the z-direction is adjusted so that the cross-section of the tube decreases as z increases. As the 
atoms move forward in the magnetic guide, the ones with a large transverse energy are evaporated. We 
then rely on elastic collisions between the remaining atoms to decrease the longitudinal velocity width and 
to increase the phase-space density of the beam. 

We investigate theoretically the efficiency of this evaporative cooling scheme by solving the Boltzmann 
equation with two different approaches: one is based on a truncated distribution ansatz for the phase-space 
distribution function, while the other one is a direct Monte-Carlo simulation of the atomic dynamics inside 
the guide. For a suitable geometry of the evaporation scheme, the results of the analytical method fit quite 
well those of the Monte-Carlo simulation. They show that, with reasonable experimental parameters, the 
cooling process can lead to a spectacular increase of the phase-space density of the atomic beam. This might 
be considered as a possible realization of a continuous atom laser ^J^. Such a coherent source of atoms 
would have fascinating applications in atom interferometry and holography, metrology and atomic clocks, 

and nonlinear atom optics 

The paper is organized as follows: in section || we describe the magnetic guide (Sec. II A), we give some 



typical parameters of the atomic beam so urce (Sec. |lIB| ), we introduce the Boltzmann equation (Sec. II C ) 
and we detail the evaporation scheme (Sec. |ll E)| ) ; in section III we use a truncated Boltzmann distribution 



to obtain an approximate analytical solution of the Boltzmann equation for the case of a ID evaporation 



scheme ; in section IV we explain the Monte-Carlo simulation and we analyze the results of the two methods 
in section 0. The paper is concluded in section ^ with a brief discussion of the coherence properties of the 
beam when the cooling is sufficiently efficient to reach quantum degeneracy, that is when the spatial density 
is of the order of A~^, where A is the local thermal wavelength. 



II. MODEL CONSIDERED 
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A. The atomic guide 



We assume that the magnetic guide consists of four parallel wires carrying the same current ±1 along the 
z direction (see fig. |l|); each wire is at a distance a from the line x ^ y = 0. The resulting magnetic field is 
B = b'{x,—y,0), with b' = 2fioI /{ira'^). We superimpose a longitudinal magnetic field Bq along the z axis 
so that the modulus of the total magnetic field can be written for x, y sufficiently small: 

B{r) = {B',+b''{x'+y')f'o.B, + ^{x'+y') . (f) 

A magnetic moment ^ which is prepared in the direction opposite to the local B will therefore experience a 
trapping harmonic potential in the transverse directions. This potential is necessarily isotropic in the x — y 
plane because of Maxwell equations for magnetostatics. The oscillation frequency in this plane is given by 

= (M6'V(mi?o))'/' (2) 

where m is the atomic mass. 

Typical experimental values are a = 4 mm, / = 500 A, and Bq = 1 mT. For rubidium atoms with /i equal 
to the Bohr magneton, this leads to an oscillation frequency ^ 2tt y. 1000 Hz. The quadratic expansion 
leading to (|^) is valid for transverse temperatures up to 500 ^K. 



B. Parameters of the atomic beam source 



In our model the atoms - for the envisaged experiment we consider rubidium atoms - enter the magnetic 
guide in the plane z = with a Gaussian velocity and space distribution with the same initial temperature Tq 
for transverse and longitudinal degrees of freedom. The longitudinal velocity distribution is centered around 
a non zero value vq. The ratio uo/Awq , where Adq = ^ ksT /m is the r.m.s. of the Gaussian velocity 
distribution, should be larger than 1 to ensure that only a very small fraction of atoms is initially moving 
"upstream" , that is with a negative longitudinal velocity. We choose in the following 

vo/Avo = 3. (3) 

For an initial temperature To =400 /iK, this gives vq = 60 cm/s. The r.m.s. of the spatial transverse 
distribution is R±q = Avo/U±. 

We further assume an initial on-axis density no at 2: = y = of the atomic beam of 8 x 10^"'^ cm^"^ which 
corresponds to a flux $0 = 2TTR^QnQVo ~ 3 x 10^ s^^ The initial stage of the evaporative cooling within 
the guide can be described by classical dynamics as no Aq w 7 x 10^'' ^ 1, where Aq = h/y/2nmkBTo is the 
initial thermal wavelength of the gas. 



C. Colhsional dynamics inside the guide 

The colhsional dynamics inside the guide can be theoretically described by the Boltzmann equation, which 
gives the time evolution of the atomic phase-space density /(r,p,t): 



dt m 



+ ^ • ^r/ - Vrf/ • Vp/ = /coU [/] , (4) 



where U{x,y,z) — Ux{x) + Uy{y) with Ux{x) = ^mil^x^ and Uy{y) = ^mfl'^y'^. The distribution / 
is normalized so that its integral over momentum gives the spatial density. Considering only elastic and 
isotropic collisions between guided atoms, we write the colhsional integral as: 

Ic,n[f] =-J^ / "i'PJ / •'"Pi j <i'Pi (/(r,Ps)/(r,P4)-/(r,p)/(r,P2)) 

.(p + P.-P3-p,)<|l + |;-ii-ii). (5) 
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Here we assumed that the cross-section a is independent of atomic momentum, a vahd assumption for 
alkaU atoms (in the absence of a zero energy resonance) if the temperature is low enough so that collisions 
essentially occur in the s-wave regime. For instance for polarized rubidium atoms, the region of s-wave 
scattering extends up to 400 /zK, corresponding to a r.m.s. velocity Awg = 20 cm/s. Above this value, 
d-wave scattering is not negligible, and it may significantly modify the results of this paper. We recall that 
p-wave scattering (more generally any odd wave scattering) is forbidden at any energy for polarized bosons. 

The collision rate of the atomic beam source 7coii is related to the on-axis density no by 7coii = 
(2/y^) no CT Awq. For the s-wave collisional cross-section of rubidium (cr = 7.6 x 10""'^^ m^) we get: 

7cou/r!± = 0.02 , (6) 

which means that an atom performs on average several transverse oscillations between two collisions, that 
is we are in the coUisionless regime for the transverse degrees of freedom. 



D. Evaporation scheme 

Evaporation along the beam eliminates atoms outside a section in the transverse xy plane. In order to 
optimize the efhciency of this evaporation, we assume that the size of this section varies with z in a controlled 
way (forced evaporation). In practice this can be achieved using a radio- frequency field at a fixed frequency 
v-c.f., and a spatially varying bias field and gradient 6', producing thus a spatially varying bottom of the 
magnetic well while keeping a constant il^ cx h' j ^/Bq [0. The stability constraints concern mostly the bias 
field Bq as it controls directly the effective local trapping potential depth A[/o — {ht^i.i./g) — M^o where g 



is the Lande factor of the atomic level. As we shall see in section § V A a typical value for AUq at the point 
where quantum degeneracy is reached is 1 fiK. This value is comparable to the potential depth typically 
used in Bose-Einstein condensation experiments, so that we do not anticipate any particular difficulty with 
the control of Bq ■ Another possibility to achieve evaporation is to place an absorbing material at a controlled 
distance Aevap(-2) from the center of the guide. 

Evaporation is taken into account in the model by putting /(r,p,<) = if the phase-space cell {r, p} 
is outside the domain where atoms are trapped. For the case of a radio-frequency evaporation, which is 
cylindrically symmetric as long as gravity is negligible, the evaporation criterion is: 

x'+y'>Al,^{z) (7) 

where Acvap(^) is determined by the radio- frequency field: g^i^B^ -\- V'^K^^^^Y^'^ = hiy^f. This corresponds to 
a 2D evaporation scheme. We also will consider in this paper a ID scheme with the criterion 

> ALp(2) ■ (8) 

As we will see below, this ID scheme allows an approximate analytical treatment. We also define the cut-off 
energy: 

ee(z) = ^mnlAl^^piz) . (9) 



III. ID HYDRODYNAMIC APPROACH 
A. The truncated Boltzmann distribution ansatz 

In this section we restrict the discussion to the ID scheme of Eq. (^). We now assume that the flow of the 
gas is in the hydrodynamic regime along the longitudinal direction, that is the macroscopic quantities - such 
as density or mean velocity - vary slowly with z on a scale given by the mean free path d{z) = Av {z)/"f coii{z). 
The assumption of a hydrodynamic regime implies a local thermodynamic equilibrium characterized by a 
local temperature T{z) at abscissa z. 

We also assume that the mean free path is much larger than the transverse extension R±{z) of the beam, 
or equivalently 7coii(2) <C fl±. As a consequence, in the ID evaporation scheme, if an atom emerges from a 
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collision with a kinetic+potential energy along the x-axis higher than the cut-ofF ec{z), it will fulfil after the 
further oscillation - and before it undergoes another collision - the condition (j^) and it will be evaporated. 
Therefore it is equivalent to formulate the evaporation criterion either in terms of the coordinate x or in 
terms of exix,px) = Ux{x) + Px/{2m). Consequently we replace Eq.(||) by 

Imnlx^ + ^ > e,iz) . (10) 
2 Zm 

The assumption of a local thermodynamic equilibrium of the gas at each abscissa z around a mean velocity 
v{z) = p{z)/m suggests the following ansatz for the classical phase-space distribution \ ^a\ : 

/(r,p) = fo{z) e-^(-)(^-+^=') e-'3(-)(p-p(-))V2™ Y{e,{z) - e,(x,p,)). (11) 

Here (3{z) = l/(kBT{z)) and Y is the Heaviside step function. The truncation takes into account the 
evaporation: f — for an atom whose energy in the transverse direction x exceeds the local depth ec(z). 
The on-axis phase-space density foiz) is calculated by the normalization condition J dxdy J (Pp f = piin(z); 
here pnniz) is the linear density of the gas, that is the number of particles per unit length in the guide. A 
straightforward calculation leads to 

f (,) ^ 1 1 Plin(^) MO] 

' l-e-'7(^) (27r)5/2 {mAv{z))3Rliz) ^ ' 

where we have introduced the thermal velocity Av{z) — y^kBT{z)/m, the thermal transverse size of the 
beam R±{z) = ^JkBT{z)/mVt\, and the quantity 

ry(z) = e,(z)/fcBT(z), (13) 

which is a crucial parameter to control the efficiency of evaporation. 

This reformulation of the evaporation criterion in the energy domain is difficult to extend to the 2D 
scheme. It would require the further assumption of ergodicity of the transverse motion, which is not correct 
for the present axi- symmetric potential since the knowledge of the total transverse energy €!_ — Cx + is 
not a sufficient criterion to determine whether an atomic trajectory will be evaporated. For instance atoms 
moving along transverse linear trajectories are evaporated when e± > £ciz), while atoms with transverse 
circular trajectories remain trapped as long as e± is below 2ec{z). This 2D case will be treated using a 
Monte-Carlo technique in section VB p^. 



B. Hydrodynamic equations 



We now derive a closed set of partial differential equations for the three macroscopic quantities T{z,t), 
Piiniz^t), and p{z,t). We note that /(r,p) in the transverse direction depends only on e^, and ey, quantities 
preserved by the Hamiltonian evolution. Therefore Eq.(H) reduces to 



df Pzdf _ 

^ H 'E~ — -'coU IJ\ 

at m oz 



(14) 



We now proceed in a way analogous to the derivation of the standard hydrodynamic equations from the 



Boltzmann equation: multiplying Eq.(|14D by 1, pz, and 
the five variables x,y,px, Py,Pz gives: 



+ Pi/"^^! respectively, and integrating over 



d d 
dlP'^-^d-z 



■ Pii 



d_ 

dz 



p- 



ksT H pii 



d_ 
di 



d_ 

dz 



= - r piin 
p- 



= - [T ^+T,kBT] p,, 
2m 



(15) 
(16) 

(17) 



where 
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^ I dxdy J d^p e,/(r,p) = fc^T (^1 - 



(18) 

Equations ( |15| , |16| , |17D are the equations of conservation for the number, the momentum, and the energy of 
the particles. On the right-hand side we have source terms due to evaporation: T is the evaporation rate at 
abscissa z, and kgT is the rate for the decrease of energy in the local reference frame moving at velocity 
v{z). 

The ID evaporation model allows to derive explicit expressions for T and Fg. We obtain for the loss 
rate of particles: 

riz) = apUz)^^^e-^^'^Sirjiz)) (19) 

The analytical expression of the positive dimensionless coefficient S{ri) is given in the appendix. As shown 
in fig. 1^, it is a nearly constant quantity, of the order of 0.075, when rj varies between 2 and 10. Eq. (p^) 
shows that the decay rate F of the linear density, given in eq.(^5|), is proportional to the coUisional cross 
section a and to the local atomic density, as expected for binary collisions. 
For the loss rate of energy, we find: 

F,(z) = F(z) (viz) + l + Siviz))^ ■ (20) 

The expression of the positive dimensionless coefficient S{ri) is also derived in the appendix, and plotted in 
fig ||. It increases from 0.43 to 0.66 when r] varies from 2 to 10. 

On the left-hand sides of equations ( p^ , p^ , pT| ) we have neglected the terms arising from the z-dependence 
of the cut-off ec{z). These terms would account for spilling, that is the loss of particles even in absence 
of collisions due to the lowering of the cut-ofF ec{z) along the beam. Neglection of spilling is valid for 
7^ 3> 1, which is well verified for the optimal evaporation with realistic initial parameters; for instance, for 
the simulations presented in section IV, we have chosen 77 ~ 5 ||lj]. For consistency we also replace in the 



following ex and /o by their values for 77 — > 00: 



(27r)5/2 {mAv{z)fRl{z) 
We finally obtain from Eqs. (|ll,0,0) in the hmit of large ry: 



Ix = keT (21) 



d p d 
ot moz I oz 



-Fpiin (23) 



It + ^Ir) P + —^AkBT pun] = (24) 

kBT{r~\v,) (25) 
5 



dt mdz J piin dz 

d P d \ , ^ "2 , ^ d [ p 
■7r^ + -j^]kBT+-kBT— ^ 
ot moz J 5 oz Im 



C. Stationary regime: z-dependence of phase space density 

We are mainly interested in a stationary regime for which we obtain a set of non-linear equations which 
can be solved by standard numerical methods: 

(Plin) - -^Plin F + F, 4^ (26) 



z P \ ^ 

^ (-) = '-f- (2^) 

oz \m/ h 

d , rnksT (2 2 fcijTX , , 
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where we have put 



5p2 7 
E=-'---kBT. 
2 m 2 



From the above we see that if p^/2m is bigger than (7/10)fcsT {i.e. E > 0), which is indeed the parameter 
range studied in this article, the mean velocity increases as a function of z. An interpretation of this result 



will be given in section V A 



A figure of merit for our scheme is the degree of increase of the phase-space density along the beam. As 
a measure of this increase we take the on-axis phase-space density fo{z) oc piin/C^s?")^/^: 



d_ 

dz 



In 



(fcsr)5/2 



= (7y + 5(r?)-2) ^r. (29) 



Since the quantity S{t]) is positive (see fig. ||), we conclude that the phase-space density increases when 

■n>2. 

IV. MONTE-CARLO SIMULATIONS 

In this section we compare results from the approximate analytical ansatz of the previous section with 
results of a Monte-Carlo simulation. The fact that the Monte-Carlo simulation requires a long computing 
time restricts the parameter space which can be explored. In particular, as explicited below, the length 
of the guide influences the memory requirement of the Monte-Carlo simulation. The total length L of the 
system, expressed in units of the mean free path do = \/t^ / {'^.n^a) is chosen in the following as 

i/do = 2500. (30) 



For the parameters given in section II B this corresponds to a length of 3.7 meters; the average time T = L/vq 
for an atom to travel from the entrance to the exit of the guide is T = 4 x 10^ (in absence of collisions 
and evaporation); this time corresponds to 830 ^^^n- 

A. Principle of the Monte-Carlo method 

This method has originally been introduced in the context of molecular dynamics [p^ . For the case of 
dilute gases, it relies on the idea that one can separate the description of the collision from that of the motion, 
allowing a simulation of the dynamics on a time scale shorter than the mean time between two collisions. 

In essence the approach consists in solving (^ numerically by evolving macro-atoms, each of which rep- 
resenting i real atoms. The macro-atoms evolve in the same potential as the real atoms, and they have the 
same initial velocity and position distributions. Their coUisional cross-section is and their initial spatial 
density is no/£, so that the collision rate 7coii and the two dimensionless parameters "Do/Auq and L/d^ are 
not changed. This approach is valid since the Boltzmann equation is invariant under the scaling a — > £(t, 
f ^ f / L Using the symmetries of the problem we restrict the evolution to the first quarter (a; > 0, ?/ > 0) 
of the entrance plane z = reducing the memory requirement by a factor 4. We evolve the macro-atoms in 
this first quarter with refiecting walls at the planes x = and y = 0. 

We inject on average 84 macro-atoms (21 in the first quarter) every Vl'^ . We take I — 5600 to match the 



flux of the atomic beam source presented in section II B, In absence of evaporation, Af = 21x4xl0'* = 



8.4 X 10'^ macro-atoms are present on average at a given time since, as stated above, it takes on average 
4 X 10"' il^^ for a particle to travel along the guide. 

Binary elastic collisions are taken into account using a boxing technique |l^-|l^ . We introduce in position 
space a lattice with a unit cell volume 5V , chosen such that the average occupation pocc of any cell is much 
smaller than 1. Collisions occur only between two macro-atoms occupying the same cell, and the time step 
5t is adjusted in such a way that the probability pcoU of a coUisional event during 5t is also much smaller 
than 1. We choose typically pocc ~ Pcoii ^ 10%. 

Evaporation is implemented in the simulation by eliminating the macro-atoms whose coordinates fulfil 
the chosen condition of evaporation. In particular, we have treated both the case of evaporation with a 
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ID position cut (Eq. (g)) in order to compare with the results of the analytical ansatz, and the case of 
evaporation with a 2D position cut (Eq. (|^)). 

We let the simulation evolve until a steady-state is reached. The corresponding time is ^ 2L/vo — 
8 X 10^ ri]]^. This allows to obtain the average energies along each axis and the linear density at a given 
location z. From these quantities we can predict the decrease in temperature, the loss of particles and in 
consequence the phase-space density increase. 



B. Comparison with the hydrodynamic approach for ID evaporation 

The results of the Monte-Carlo and of the hydrodynamic approaches for ID evaporation are plotted in 
Figs. 1^,^ and giving the variations with z of the mean velocity v{z), the flux ^{z), and the phase-space 
density foiz). This set of data has been obtained by choosing an evaporation barrier edz) approximately 5 
times larger than the mean local transverse energy, that is we kept ry ~ 5; as shown below, this value for rj 
leads to a gain of 7 orders of magnitude for the phase space density with a minimal length requirement. 

Using the Monte-Carlo simulation, we have found that the three components of the velocities have a nearly 
Gaussian distribution, with dispersions equal to within a few percent. This is consistent with the hypothesis 
of local thermodynamic equilibrium at the basis of the hydrodynamic approach. 

The results shown in figs. I, I and I show an excellent agreement between the hydrodynamic approach 
and the Monte-Carlo simulation. This allows one to make all optimization and design procedures for the 
choice of initial parameters and spatial variations of the cut using the approximate hydrodynamic treatment 
while keeping the Monte-Carlo simulation - which requires several days of computation on a work station - 
for final checks. 



V. DISCUSSION OF THE RESULTS 



A. ID evaporation 

Using the equations ( p6[^ ) in the regime of evaporation with a constant ry, we have first determined the 
optimal choice for rj. As shown in fig. ^, the shortest distance providing a lO'' gain in phase space density 
- required in order to reach quantum degeneracy for the initial conditions considered in §2- is obtained for 
T] ^ 5. If smaller phase space gain are needed, smaller values of ry would be more appropriate since they 
would lead to shorter cooling lengths. 

We now consider fig.^, obtained for rj — 5, which shows a slight increase of the mean velocity v{z), as the 
atoms progress within the guide. This increase can be understood from kinematic arguments: the cooling 
along z can be seen very crudely as a process in which a fast and a slow particle collide, one of them being 
eliminated. Now, at a given location, the fast particles are renewed with a larger rate than the slow ones, 
because it takes less time for them to go from the entrance plane to the considered location. Therefore, the 
mean velocity in a location z > 0, where particles have already undergone in average several collisions, is 
larger than in z = 0. This acceleration effect becomes negligible as soon as kBT{z) <C p^/{2m) since the 
beam is then quasi mono-kinetic. 

The figs. show the flux and the phase-space density as a function of z for 77 = 5. A very significant 
increase of phase-space density by a factor 500 proves the efficiency of the evaporative cooling. This increase 
is accompanied with a reduction of the flux by a factor 5.5. The phase-space density and the flux both vary 
quasi-exponentially with the position z. 

To get a better understanding of the variations with z of these quantities, we now consider Eqs. ( ^6| - ^ ) 
in the limiting case where v is constant and fc^T <C E. Assuming also a constant rj, the solutions of these 
equations are: 

Ain(^) = Plin(O) (l - ^\ T{Z) = r(0) i^-^X (31) 



with 



V Se" 5 2(7? + S{t^) - 1) 

dn -; ; = (Xc — = Pc = (32) 

Az;o 7r3/2S'(7y)(7y-hS'(r/)-6) + 5(r/) - 6 r? + ^(r/) - 6 
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For the parameters of figs.^,|| and |[ where = 5, we find ac — —12, /3c — —21.9 and Zc — 
—4300 dov/Avo ~ —13 000 do, so that the relevant range of lengths z are much smaller than \zc\- Con- 
sequently we can approximate the previous results with: 

Piin(z) ~ piin(O) exp(-Q;cz/zc) r(z) ~ r(0) exp(-/3c2/zc) (33) 

This leads to exponential variations exp(z/zp.s.) and exp(— z/zgux) of the phase space density and the flux, 
with: 

This explains why the results in figs. exhibit a quasi-exponential variation with z. For 77 = 5, we get 
Zp.s. ~ 100 dow/Awo and Zflux — 360 dow/AwQ. The factor 3.6 between these two lengths suggests that a 
phase space gain of 7 orders of magnitude can be achieved with a flux reduction by less than 2 orders of 
magnitude. This is confirmed by the numerical solution of the hydrodynamical equations Eq. ( p6[p8| ) : the 
required length is 7600 do, on the order of 11 meters for the parameters of section it is accompanied by 
a decrease of temperature by a factor 4000 and of flux by a factor 90. 

It is worth noting that the collision rate 7coii(z) does not vary much for z <C |zc|- Indeed, we have: 

7coii(^:) = 7cou(0) 1 = 7coii(0) — ■ 

The sign of Zc is the same as the sign of 77-I- S'(?7) — 6, which vanishes for 77 ~ 5.4. Therefore, the collision rate 
7coii increases with z if 77 > 5.4 and decreases if 77 < 5.4. For the particular value 77 — 5.4, the collision rate 
is constant, as well as the loss rates for particles r(z) and energy Tf{z); the quantities p\in{z) and ksTlz) 
then have an exponential variation for any z, since Zc,ac, Pc — ^ 00. 



B. 2D versus ID evaporation scheme 

The experimental setup based on radio-frequency evaporation and discussed in section ||, corresponds to 
a 2D evaporation scheme while the analytical treatment of section III is based on a ID evaporation. One 
does not expect the two situations to be equivalent. More precisely, since the x-y degrees of freedom are 
not mixed by the coUisionless motion in the axi-symmetric potential, the ID evaporation is expected to be 
less efficient than a 2D evaporation scheme: a particle may emerge from an elastic collision with a large 
transverse kinetic energy along the y-axis, without being evaporated in the ID evaporative scheme. 

We have checked with our Monte-Carlo simulation that 2D evaporation is more efficient indeed than ID 
evaporation. We have first used a crude model of 2D evaporation, with the truncation of Eq.(pT|) replaced 
by Y{ec{z) — etot + kBT(z)/2), where etot = e± + ipz — p{z))'^ /2'm is the total energy in the frame moving 
at velocity p{z)/m. This truncation assumes in particular quasi-ergodicity in the x — y plane. Within this 
model we perform the same optimization as in fig. ^; to achieve a gain in phase space density by 7 orders of 
magnitude the optimal 77 is now equal to 6. We have then run the Monte-Carlo simulation for this value of rj. 
The gain of a factor 500 in phase space density, which was obtained for a length 2500 do for ID evaporation 
(see fig. H), is reached now for a length 900 do; the relative variation of the flux (reduction by a factor ~ 6) 
is similar. The variations with z of the phase space density and the flux are quasi-exponential. Note that 
the predictions of the crude model are in amazingly good agreement with the Monte-Carlo results. 

To conclude the length of the guide needed to achieve a specifled gain of the phase-space density is a factor 
of ^ 3 smaller for 2D than for ID evaporation. E. g. the 11 meters of evaporation length necessary in ID 
for the experimental conditions of §0 are now reduced to 4 meters. 



VI. PERSPECTIVES 



For a typical experimental source of cold atoms, the length needed for a phase-space increase of seven 
orders of magnitude ~ which should bring the system close to the degeneracy point - is of the order of a few 
meters. The possibility of quantum degeneracy in such a system raises interesting questions. As it is well 
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knownM], there is no Bose- Einstein condensation in the thermodynamic hmit in a ID geometry, obtained 
here ]20[ in setting A^, L — > oo, while keeping a constant linear density pn^ = N/L, a constant temperature 
T and a constant transverse oscillation frequency il± . Therefore we do not expect a macroscopic occupation 
of a single quantum state of the longitudinal motion. 

In order to get more insight in the output of this system, we assume that the transverse extension of 
the thermal component of the beam R± = ^ksT / {mVi\) is much larger than the thermal wavelength A = 
h/^fhrmk^ (that is /cbT ^ h^±). In this case, elastic collisions ensure that thermodynamic equilibrium is 
reached in the frame moving with the mean velocity of the gas. In this frame, we describe the properties of 
the system using as an approximation the grand canonical Bose-Einstein distribution for an ideal gas. The 
gas being at temperature T, we obtain the maximal linear density that can be put in states corresponding 
to an excited transverse motion 

with C(5/2) — 1.34. When one increases the linear density above this critical value (which corresponds 
to a spatial density on axis larger than n^'^'> = C(3/2)A~'^), the transverse degrees of freedom undergo a 
Bose-Einstein condensation J2^ . that is atoms start accumulating in states corresponding to the ground 
transverse motion (see fig. l7|). By convention, we set the energy of this ground state to zero. Since the 
chemical potential n then satisfies <C h^±_ <^ ksT, the longitudinal momentum distribution for these 
atoms can be approximated as: 

"^^^ " h exp[(p2/(2m)-^)/fcBT]-l ~ ^ + pi ^^""^ 

that is a Lorentzian distribution of half width pc = [—2mfj,]^^'^. Such a distribution leads to a spatial 
correlation length of the gas along z given by = fi/Pc- By integrating n{p) over p, we can relate fi and 

(c) 

therefore to the linear density of atoms in the transverse ground state pn^ — p\^^ : 

, 1 , , , h n — n(^) , , 

^c = ^(.n„-.,t^)=2C(3/2);^^^ (37) 

where n is the 3D on-axis density. This correlation length can be much larger than the thermal wavelength. 

(c) 

Thus, for piin > Piin , the output of the system can be viewed as propagating independent trains of matter 

waves, each of which having a length of the order of and containing on average A(piin — /(27r) 
atoms. 

To conclude, we have presented in this paper the principles of the evaporative cooling of an atomic beam. 
In the classical regime, where the mean interparticle distance is much larger than the thermal wavelength 
A, we have evaluated the characteristic lengths for the flux and the phase space density variations. We have 
also outlined briefly the coherence properties of the beam once quantum degeneracy is reached. A more 
detailed characterization of these coherence properties, including the effect of interactions between particles, 
will be the subject of a future work. 
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VII. APPENDIX 



In this appendix we derive explicit expressions for the loss rates for particles F and for energy in Eqs. 
( |l5| , p^pT| ) . Integrating over x, y, and p to derive (^5|), we obtain: 

F(z)piin(z) = JdxdyJ (fp j dPp2 j (fpz j d^Pi (/(r, p) /(r, pz) - /(r, pa) /(r, P4)) 

^(P + P2 -P3 -P4) <5 ffl + fi- - fi- - . (38) 
\2m 2m 2m 2m J 

In this integral, the phase space cell (r, p) is the variable of integration of the left hand side of (|^) and it is 
always in the trappable domain defined by: 

|i- + ImQlx' < e,(z) (39) 
2m 2 

For the first part of the integral, representing the collision p + P2 ^ P3 + P4 and involving /(r, p) /(r, P2), 
the cell (r, P2) is also in the trappable domain, while the two cells (r, pa) and (r, P4) may be either in or out 
of the trappable domain. For the second part of the integral, representing the collision P3 + P4 — > p + P2 
and involving /(r, pa) /(r, P4), the two cells (r,p3) and (r,p4) are in the trappable domain, while the cell 
(r,p2) may be either in or out of the trappable domain. 

We now rearrange the second part of this integral by exchanging the role of (p,P2) and (p3,P4). After 
cancellation of various terms, we are left with: 

F(z) piin{z) = J dxdy J d^pi J d^p2 J d^ps J d^p^ /(r, pi) /(r, P2) 

.(p,+P.-P3-p,)*(ii + ii-|;-i|). (40) 

where (r, Pi) and (r,p2) are in the trappable domain, (r,p3) is out of the trappable domain, and (r,p4) is 
either in or out of the trappable domain. 
We now put 

p ^ P1+P2 ^ Pi P2 p, ^ P3+P4 , ^ P3 P4 

2^2 2^2' ^ ' 

The (5-distributions entering into ( ^0|) impose P = P' and |q| = |q'|. We perform some rearrangements, 
using as integration variables for the vectors q and q' the coordinates QxtQx along x, the moduli q,q' and 
the azimuthal angles around x-axis. We integrate over those angles and we split the integration domain into 
a part where only particle 3 escapes and a part where both 3 and 4 escape. We then obtain: 

/Xc pQ roc />oc />oc 

dx dPx dy dPy dP, e~"^U.+U,^/k.T ^-Pyrnk.T 

-Xa Jo J — OC J —00 J — OO 

dqx I dqq e'" /"'^-^ / dg^ 



j-ca r pQ+Pj, pq 

/ dq q e-9 /"'^'^^ / dq'^ + 2 dq', 

JQ+P^ JQ-P. JQ+P^ 



Q+P^ rq 
dq'^ + 2 
Q-P^ Jq+P=o 



(42) 



where P — |P| and q — |q|. The length Xc{z) is the positive solution of m{l\x1{z) — 2ec{z), and Q{x,z) is 
the local escape momentum given by /2m = ec{z) — Ux{x). The expression (^) is the direct transcription 
of Eq. (26) of [|l3|, to the situation considered in the present paper. 

The integrations over y,Py,Pz are immediate, as well as for qx- After an integration by part over q, we 
get: 



dw e-'"^ / du{a- u) e""' (erfc(a - u) + erfc(a + u)) (43) 
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where we have put w = x/R±, a = Q/\/mkT = y^2ri — w"^, u — Pxj^ mkT and: 

2 



erfc(u) 



dv e"" 



Using finally the relation (M/m to express /o in terms of the linear density pun, we arrive to (|19|) with: 



27r2 (l-e-'?)2 J_ 



V2v 



dw e ^ / du {a — u) e (erfc(a — w) + erfc(a + u)). 



(44) 



The calculation of proceeds along the same line. We multiply (Q) by tx + f-y +p1/2m and integrate 
over x, y, and p. Using the same rearrangement as above, we get: 

/ dPx / / dPy / dP. e-2(^^+^'')/'=-^ e-PV-^-..T 

•'0 J — oo J —oo J — oo 



dq q e 



— g /7nkBT 



E3 dq'^ 



dq q e 



^q^/mkgT 



Q+P^ 



fQ+P^ f-q 11 

/ E3dq',,+ / iEs + E,)dq'^ \ 
JQ-P. JQ+P^ J J 



(45) 



where: 



^4 = 1^ + C/X + C^, = ^^'^ +Ux + Uy. 

2m 2m 



The expression ( p5[ ) is also the direct transposition of Eq. (26) of |13|. After integration and changes of 
variables similar to the ones given above for the calculation of 5* (77), we reach (pOl) where 



47r2 5(?7) (1 - 6-^)2 J_ 



dw e ^ / du {a — u) e ^ 



+ ^__^e~("+")' - (a^ - u2)(crfc(a - m) + crfc(a + u)) 



(46) 



with a = y^2r/ — w"^ as in ( p4| ) 
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FIG. 1. An atomic beam propagates in a transverse magnetic guide. Evaporation eliminates particles whose 
transverse coordinate exceeds some adjustable value Aevap(z). The emerging beam is colder and it has a larger 
phase-space density than the input beam. 

FIG. 2. Variations with of the dimensionless parameters S multiplied by 10 for clarity (full line) and S (dashed 
line). The analytical expressions for these parameters are derived in the appendix. 

FIG. 3. Longitudinal velocity as a function of the position z measured in units of the initial mean free path 
do = / i^rioa) . The unit for velocity is the initial velocity spread Avq- The continuous line corresponds to the 
numerical solution of the ID hydrodynamic equation for rj — 5. The markers indicate the results of the Monte-Carlo 
simulation. 

FIG. 4. Flux as a function of the position z measured in units of the initial mean free path do ~ \/TT/{2nQa). 
The unit of flux is the initial value ai z = 0. The continuous line corresponds to the numerical solution of the ID 
hydrodynamic equation for 77 = 5. The markers indicate the results of the Monte-Carlo simulation. 

FIG. 5. Gain in phase-space density as a function of the position z measured in units of the initial mean free path 
do ~ \/n/{2nQa). The continuous line corresponds to the numerical solution of the ID hydrodynamic equation with 
■q = 5. The markers indicate the results of the Monte-Carlo simulation. 

FIG. 6. Gain in phase space density as a function of the position z expressed in units of the initial mean free path 
do = -^^/(Swocr). These curves are the numerical solutions of Eqs. {bq-psl) in the case where rj is fixed. 



FIG. 7. For an ideal Bose gas fraction of atoms in the transverse ground state, as a function of the linear density. 
This figure has been obtained for ksT = 20 hi^±, so that p^^^ ~ 536 A^^. 
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